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MATRIX METHODS FOR DETERMINING THE LONGITUDINAL-STABILITY 
DERIVATIVES OF AN AIRPLANE FROM TRANSIENT FLIGHT DATA 1 

By James J. Dokesan 


SUMMARY 


Three matrix methods are 'presented for determining the 
longitudinal-stability derivatives from transient flight data. 
One method, which requires four measurements in time-history 


form and utilizes the incremental tail load to separate the pitch- 
ing-moment derivatives C m g and C m£l , permits the computation 


.of all the longitudinal-stability derivatives. A second method 
| requires three measurements and one supplemental assumption, 

! C m . 

\namely Constant. This method gives the most infor- 

L ' ma 


| mation for the least amount of work. The third method re- 
quires two measurements and two supplemental assumptions, 


C 

namely Constant and C, H =-A C L . ( where C ms and C L} 


are the elevator-effectiveness derivatives, x t is the tail length , 
and c is the mean aerodynamic chord). An inspection of the 
results obtained for the various methods shows the scatter which 
is typical of this type of analysis of flight data. 


INTRODUCTION 


The determination of the longitudinal-stability derivatives 
from flight data is a relatively difficult task because the 
wind-tunnel technique of permitting only one variable to 
change at a time, while constraining all the rest of the 
variables, cannot always be used. It is in the analysis of 
such flight-test data that matrix techniques employing the 
equations of motion seem to be particularly useful. 

Currently, much work is being carried out on the determ- 
ination of stability derivatives directly from flight data but 
as yet this work is still in the preliminary stages. The 
matrix methods for the determination of stability derivatives 
from transient flight data that are developed herein are an 
addition to this work. The previous work done on the 
determination of longitudinal-stability derivatives is exten- 
sive, and no attempt is made to summarize it since this 
summarization has been adequately done in reference 1. 

In the present report three methods are developed and 
presented for determining the longitudinal-stability deriva- 
tives from transient flight data. In these, methods the 
expressions for some of the stability derivatives are in the 
form generally used in stability calculations. The first 
method requires the combination of four measurements in 
time-history form, two of which must be incremental elevator 
deflection and incremental tail load and the other two 


measurements can be chosen from a possible three, namely 
incremental load factor, pitching velocity, and angle of 
attack. The method demonstrates the use of the tail load 
to separate the pitching-moment derivatives On and C m& 
and to determine the downwash derivative be/ba. 

The second method, which is more restricted, requires a 
combination of three measurements (in time-history form), 
one of which must be incremental elevator deflection and the 
other two measurements can be chosen from a possible three, 
namely incremental load factor, pitching velocity, and angle 
of attack. This method also requires one supplementary 
assumption, namely C m& =\C m g, where X is a constant! 

The third method uses a combination of two measurements 
(in time-history form), one of which must be incremental 
elevator deflection and the other one may be chosen from 
incremental load factor, pitching velocity, angle of attack, 
and so forth. The method also requires two supplementary 

assumptions, namely C m .=\C m j and C ms —--~-C Li (where 

c 

C ms and C L( are the elevator-effectiveness derivatives, x t is 
the tail length, and c is the mean aerodynamic chord). By 
using a modification of the third method, it is shown that 
considerable information can be obtained from a single time 
history. 

The methods are demonstrated by applying them to flight 
data obtained from tests of a medium jet bomber, and a 
comparison of the derivatives obtained .by the various 
methods gives an indication of the accuracy which can be 
expected from data analysis by matrix techniques based on 
the longitudinal equations of motion. 

SYMBOLS 


b wing span, ft 

c mean aerodynamic chord 

Ci, C 2 constants defined in appendix E 

C L lift coefficient, L/qS 

C La rate of change of airplane lift coefficient with 

angle of attack per radian; see appendix E 
C Ls rate of change of lift coefficient with elevator 

deflection per radian; see appendix E 
C L& rate of change of lift coefficient with a per 

radian; see appendix E 

C h g rate of change of lift coefficient with pitching 

velocity per radian; see appendix E 


1 Supersedes NAOA T.\ 2902, ‘‘Matrix Methods for Determining the Longitudinal-Stability Derivatives of an Airplane From Transient Flight Data” by James J. Donegan, 1953. 
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pitching-moment coefficient of airplane, 

M/qSc 

rate of change of pitching-moment coefficient 
with angle of attack per radian; see ap- 
pendix E 

rate of change of pitching-moment coefficient 
with elevator deflection per radian; see 
appendix E 

rate of change of pitching-moment coefficient 
with pitching velocity per radian; see 
appendix E 

rate of change of pitching-moment coefficient 
with a per radian; see appendix E 
pitching-moment coefficient of horizontal tail 
surface, M t jqSfit 

acceleration due to gravity, ft/sec/sec 
airplane moment of inertia, slug-ft 2 
airplane radius of gyration about pitching 
axis, ft 
lift, lb 

airplane mass, W/g, slugs 
pitching moment of airplane 
airplane load factor 

rro 

dynamic pressure, lb/sq ft 

■wing area, sq ft 
horizontal-tail area 
time, sec 

true velocity, ft/sec 
, airplane weight, lb 

length from center of gravity of airplane to 
aerodynamic center of tail (negative for 
conventional airplanes), ft 
K e coefficients of transfer function relating 8 
and 5; see appendix E 
wing angle of attack, radians 
tail angle of attack, radians 
flight-path angle, radians 
angle of pitch, a +7 
elevator deflection, radians 
downwash angle, radians 


C„ 


C mi 


X= 

Vi 
p 
r 

Matrix notation: 


tail efficiency factor, qjq 
mass density of air, slugs/cu ft 
dummy variable of integration 


|| || rectangular matrix 

[ ] square matrix 

{ } column matrix 

||Cj| integrating matrix (see table I) 

||j?,|| \\D\\, 1 1£| | rectangular matrices defined in appendix E 
Subscripts: 

i denotes row elements in matrix 

t tail 

WB wing-body combination 

Tor sign conventions used, see figure 1 . 

A dot over a symbol denotes the first derivative with 
respect to time, and two dots over a symbol denote the 
second derivative with respect to time. 


ec.g. 




Figure 1. — Sign conventions employed. Positive directions shown. 

i 

The symbol A refers to an incremental value. Inter- 
mediate variables such as Ap, A?, Act, A<p, and Aip and the 
constant K 10 are defined in appendix E. 

OUTLINE OF METHODS 

The three methods are based on the longitudinal equations 
of motion for horizontal flight and use matrix methods to 
analyze time histories of measured quantities. The equa- 
tions of motion used in each of these methods are expressed 
in the form 

^An=C La Aa+C%8+C Lit a+C Li A5 ( 1 ) 

-^=9~C ma Aa-\-C mil a-\-C m ^-\-C ms AS ( 2 ) 

These equations apply to a rigid airplane and are based 
on the usual assumptions of linearity, small angles, and no 
loss in aii-speed during the maneuver. The equations are 
further restricted to the range in which the variation of the 
derivatives is linear and also to conventional wing-tail 
configurations in which the major contribution to damping 
in pitch is due to the horizontal tail. All the variables 
are given in incremental form measured from a steady - 
flight trim condition. 

As indicated in reference 2, the four values Aa, a, and 
A 8 in equations (1) and (2) are linearly dependent; therefore, 
if four simultaneous^equationsjare formed to determine 
either the force or moment derivatives, they cannot bo 
solved uniquely for the unknowns. 

For purposes of analysis the moment equation (2) is 
integrated once and expressed in the form 

e=C ma £ Aa dt+C m .A a +C mi A8+C ms f‘AB it (3) 

This form permits the use of numerical integrating methods 
that are more desirable than numerical differentiating 
schemes when applied to flight data. Integrations of the 
variables are performed by use of an integrating matrix 
|| Cji derived in reference 3 and given in table I herein. For 

instance, a time history of A ndt may be obtained from 
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TABLE I.— INTEGRATING MATRIX \\C\\ 
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a time history of An as follows: 

| j* An dt | = He'll { An t } (4) 

The integrating matrix ||<7|| given in table I may be used 
for any time interval At; most of the computations of this 
report are based on a time interval of A#=0.1 second. This 
interval may be too large in some cases, and, if greater 
accuracy is desired, a shorter time interval may be chosen. 

The essential differences in each of the methods are in 
the number of quantities to be measured. Method A 
requires four basic measurements in time-history form to 
determine all the derivatives. Method B requires three 
measurements and one supplemental assumption, namely 
:C m& ='hC m g. Method C requires two measurements and two 

supplemental assump tions, namely C ma = \C m g and C ms = 

c 

All measurements of flight data used are time histories of in- 
cremental values measured from a trimmed level-flight initial 
position. The development of the equations for each method 
is covered in appendixes A to C; in the body of the report 
the methods are outlined by stating the pertinent equations 
in the order of computation. Since these computations make 
extensive use of least-squares procedures and are greatly faci- 
litated by the use of matrix algebra, most of the equations 
are given in matrix form. 


or 

{ a t } = { 6 t } —y { An t } 


(2) Calculate time histories of Act, Ad, 




A Sdt, 


and j J A8 dr dt by using the integrating matrix ||C|| and 
the time histories of a, 8, and AS; for example, 


{Aa<} = IICH {«<} (6) 

(3) Determine C La and C L} by least squares from the 
relation 

C La Aa+C H A8^An (7) 

or 

(4) Compute the coefficients K\, K 2 , K s , and Kt of the 

transfer function relating pitching velocity and elevator 
deflection by the use of the method of reference 4 and the 
equation ; 


. KiA9+K 2 f Addt— K s f A 8dt—K 6 ( f ASdrdt=—8 
Jo Jo Jo Jo 


( 8 ) 


METHOD A 

Of the four basic measurements required with method A, 
two must be incremental elevator angle and incremental 
tail load and two other measurements can be chosen from 
a possible three, namely incremental load factor^, pitching 
velocity, and angle of attack. In this report, incremental 
oad factor and pitching velocity are used. 

The procedure of computation with method A (see ap- 
pendix A for development) is as follows: 

(1) Compute a time history of rate of change of angle of 
attack a from 


a=9—y An 


(5) 


where the measured values of pitching velocity and elevator 
deflection are used. 

(5) Determine K w from the relation 




grSc 


K, 


(9) 


by using the results of steps (3) and (4). 

(6) Calculate time histories of the intermediate quantities 
A i p and An by using the expressions 


[t An ~xl 1 ) A “] ( l °) 

or 

{Atpi}=y^{Ani}— (V’L + f) {Aa«}., 
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and 

K w Aa-f -Kioa— ft 5 AS (11) 

or 

( t Aft, } “|t *io { *** } ~ { <* > - ft* ( AS, } 

(7) Compute C m ^ by least squares from the relation 

C m #A<p=An (12) 

or 

C m $ { A^h } = } Am } 

(8) Determine C mA from the equation 



^ m a 

(13) 

(9) Calculate the time history of the intermediate quantity 

Ao- from 

Act — B C mi Aa C m $A6 

(14) 

or 

{ A<T( }=-^<y= { 6i } C m& {Aat } C m fr{ ABi } 


(10) Compute C Ma and <7,,,, by least squares by using the 

relation 

C ma J Q Aa dt+c mi j ^ A5 dt—Ao 

(15) 

or 

||D||{g;;}-(A..) 



(11) Calculate ft(j and C La from the following definitions: 


Czi—Ont (16) 

Xt 

and 

(17) 

(12) Detennine the time history of the intermediate 
quantity Aty from the equation 

A*=pAn — Cl$ — ft^a (18) 

or 

{A^} = |§{Ani}— C Li {0i}-C Lii {ai} 

(13) Compute the refined values of C La and ft* by least 
squares by using the relation 

Ci^Aa-r CtjAS = Ai/» (19) 

or 

■nail 


These values now include the effects of the and < \ a 
terms in the force equation. 

(14) Method A can now be iterated to obtain better 
values of the derivatives by starting the process over at 
step (5) with the improved C La and C L , values from step (13) 
and following the procedure again. The iteration converges 
rapidly. 


(15) 


The derivatives ^ and 

( dCr.\ ScV p 

be 1 

yfift ft* 


are found from 


( 20 ) 


( 21 ) 


bCjL, 8 r 
55 q ,S, i{ 


( 22 ) 


This procedure shows that the derivatives Cz, a , O*, C La) 

ft J> ft‘a> ft"9’ ft**’ (~5“)/ £) a ’ ail( ^ nl! ^ 

mined by numerical operations on four time histories of 
measured flight data and through the use of the theoretical 
relationships given as equations (16), (17), (20), (21), 
and (22). 

METHOD B 


Three basic measurements are used in method B, one of 
which must be incremental elevator angle and the other 
two measurements can be chosen from a possible three, 
namely incremental load factor, pitching velocity, and angle 
of attack. In this report, incremental load factor and 
pitching velocity were used. 

In lieu of the fourth measurement, the supplemental 
assumption is made that 

C mA =\C me (23) 

If a value of X is not known in advance, a first approximation 

(see ref. 5) is X=~- Although the assumption X— ^ imposes 

a restriction on the generality of the method, it appears to 
be justified since it reduces computation time to almost 
one-half that required for method A and for the examples 
presented herein gave results which are in good agreement 
with those of method A. 

The method is outlined by merely stating the appropriate 
equations, the development of which is contained in appendix 
B. The procedure is as follows: 1 

(1) Compute the time history of Aa by using equations (5) 
and (6). 

(2) Determine a time history of the intermediate quantity ’ 
A£ from the expression 

A£=(l+X)0-yf A n (24) 

or 

{A£(}=(l+X){0j}— yf { An t } 
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(3) Calculate time histories of f A adt, f A %dt, and 

t Jo Jo 

■ J A Sdt by using the integrating matrix ||(7|| and the time 
I histories of A a, A£, and AS. 

(4) Compute G ma , C m §, and C mi by least squares from the 
relation 

c na J ' A<* dt+ Cmij'u; dt+ C ms j*A5 dt=~ 6 (25) 


or 




{M 


(5) Determine C m& from C mi by using equation (23) and 
then determine C Ljs and C L£t by using equations (16) and (17). 

(6) Calculate the time history of the intermediate quan- 
tity A^ by inserting these values of C Lh and C L g into equa- 
tion (18). 

(7) Compute the values of C La and C Ls from equation (19). 
/dCA . &Cl, 


(8) The derivatives 


V da /< 


and 


dS 


- are then deter- 


mined from equations (20) and (22) and the previously 
derived quantities. 


( METHOD C 

I Method C is an extension of the method presented in 
, reference 4. Appendix C contains the development of the 
| pertinent equations upon which method C is based. Two 
[basic measurements are used in this analysis, one of which 
'must be incremental elevator deflection and the other one 
may be chosen from incremental load factor, pitching ve- 
locity, angle of attack, and so forth. In this report incre- 
mental pitching velocity is used. 

Two supplemental assumptions are made. The first is 
the relation between O mi and C mi given in equation (23) and 
the second is 

° L ' (26) 


The procedure for method C is as follows: 

(1) Compute the stability coefficients K u K t , K b , and K 6 
as outlined in step (4) of method A. 

(2) Compute C La from the relation 

a.— (27) 

where 

Gl= W [ x S +(1 +x) (28) 

and 

(29) 


13) Determine C m by using the expression 


C m 




qS 


KJ 


qSc mV cFm(l-)-X) 


cFm(l+X) 


( 30 ) 


(4) Calculate G mi from 


Cm »~eVm{l+\) 
(5) Compute G m . from 




I mV 

n _?'ST QSc 

C L 


Ks 


O m 


x. 


(31) 


(32) 


(6) Determine C m& and C L . from equations (23) and (26) 
by using the values of G mi and G ms found in steps (4) and (5), 
respectively. Approximate equations for the stability deriv- 
atives are given in appendix C. 

In appendix D, method C is modified slightly so that 
many of the stability derivatives can be obtained from a 
single time history. This time history must be the response 
to an input elevator motion of the impulse type. This 
modified method C comes closest to the ultimate aim of 
this type of analysis, namely to determine the derivatives 
from a single time history. 

One of the important factors in obtaining reliable results 
with the methods outlined herein is the choice of a sufficiently 
small time interval At. In the computations using method 
C in this report, in one case a time interval of Ai— 0.1 
second was found to be too large to give reliable results, and 
a time interval of A<=0.05 second had to be used. 

EXAMPLES 

In order to illustrate the methods outlined in the previous 
section as well as to compare the results obtained, a number 
of examples are given in which the data used are from test 
runs of a medium jet bomber at about the same Mach 
number. Methods A and C are applied to flights 1 and 2; 
whereas all three flights are analyzed by method B. Com- 
putations are shown in the tables for flight 1, but for the 
other flights only the results are given. 

Table I contains the integrating matrix \\C\\ based on 
Simpson’s law (ref. 3) which is used in all three methods. 

The airplane characteristics and flight conditions are 
shown in table II (a) for all three flights. Although the 
geometric parameters are the same, the parameters such as 
weight, speed, Mach number, center-of-gravity position, 
and altitude vary slightly between the three runs. 

In table II (b) the coefficients of the transfer function 
which relates pitching velocity to elevator deflection de- 
fined by equation (8) and computed by the method out- 
lined in reference 4 are shown. These preliminary constants 
arc required in methods A and C and the actual computa- 
tions are shown in a subsequent table. 

Time histories of measured and derived quantities for 
flight 1 are shown in table III. The quantities in columns 
©i ©j ®) and © are measured and the other five quantities 
are derived from the measured quantities. In these tables 
more decimal places are carried in the measured quantities 
than are warranted by instrument accuracy in order to as- 
sure no loss in accuracy in rounding off. The measurements 
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of incremental tail load AL, were available only for the times 
listed, and, since these covered approximately the natural 
period of the short-pei-iod oscillations of the aircraft, the 

TABLE II.— AIRPLANE CHARACTERISTICS, FLIGHT CON- 
DITIONS, AND TRANSFER-FUNCTION COEFFICIENTS 

(a) Airplane characteristics and flight conditions 


1 Flight 1 

| Flight 2 

j Flight 3 

b, ft 

I 89 

89 

89 

c, ft I 

1 14.016 

14.016 

14. 016 
27.44 
0.062686 
258,957 
141. 61 
0. 496 
1828.66 

Center-of-gravlty position, percent M.A.C 

27.34 

0.061923 

27.32 
0.062854 
255, 276 
141.61 
0. 494 


255,865 
141.61 
0. 497 

t,», ft* 


1 1806.83 

1802. 67 

q, lb/ft 2 

171 

106 

171 

1,175 

289.2 

1 S, ft*. I 

1, 175 
289. 3 

1, 175 
289. 3 

St, ft* ! 

t .ft/sec - 

w, lb 

m <ts. j 

i,, ft 

520 
58,180 
0.289561 
— 33. 5 
0. 87 

512 
58.050 
0. 297595 
—33.5 
0.87 
0.001276 

514 
58, 880 
0.293060 
j -33.5 
0.87 
0.001281 

o, siugs/ft* :.7. j 

0. 001267 


(b) Coefficients of airplane transfer function 


1 

Flight 1 

Flight 2 

\ 

Coefficient j Probablc 
| error 

Coefficient 

Probable 

error 

1 K t 

4. 14 | 0. 14 

9.547 0.73 

-9.767 0.35 

-14.624 1.4 

4.19 
10. 329 
-10.010 
-16.526 

0.18 

0.70 1 

0.42 

1.3 | 

Kt 

. AY 

I A'a - 



data were considered sufficient. More of the time histories 
of the other variables were available and were used. 

Method A. — The principal computations illustrating 
method A are presented in table IV; some of the intermed- 
iate steps outlined in method A are simple computations and 
are therefore not included in this table. Table IV fa) is 
obtained by applying equation (7) to the data given in table 

III and illustrates step (3) of method A. ’ 

In table IV (b), the computations illustrating the deter- 
mination of C mi and C„ la by steps (7) and (8) of method A ' 
are shown. Two of the columns are taken from table III ! 
and the equations upon which the computations are based ■ 
are (12) and (13). i 

Table IV (e) illustrates the computation of C v , a and C mi ' 
by step (10) of method A. Two of the columns are ob- 
tained by operating on columns © and ® of table III with 
the integrating matrix [jfTjl'Iiven in table I, and the other 
column is taken directly from table III. The computation I 
is based on equation (15). | 

The refined values of C La and C Ls are determined in table 

IV (d) by step (13). Two of the columns are taken directly 

from table III and the other column is derived by use of 
equation (18). i 

Final results obtained with method A for the data of) 

flights 1 and 2 after three iterations arc shown in table 
IV (e). I 


TABLE III.— TIME HISTORIES OF MEASURED AND DERIVED QUANTITIES FOR FLIGHT 1 




Measured 




Derived 


i 

,© 

© 

© 

© 

© 

© 

© 

® 

® 

© 

t 

AS 

AA* 

An 

e 

Ace 


A* 

&p. 

Air 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o ; 

.1 


600 

0 


0 

M^EIHHB^B 



o ! 

.2 

.055812 

3773 

.03792 

■EWiTr^™ 


-.030886 



BHu 

.3 

.072880 

3488 


-.074064 

-. 005175 




-.007957 : 

.4 

.074625 

1644 

-.35392 





-. 022$67 

-.014935 

.5 

. 071414 

25 

-.62568 

-. 149419 



. 328389 

-.046165 


.6 

. 070698 

-1423 

-.92272 

-. 165781 


-.220103 

. 474922 


025122 

.7 

.007923 

-2529 

-1. 16920 


-.046216 

-. 215050 

. 010391 

bmH 

-.028773 | 

.8 


-3459 

-1. 40936 

-. 163628 

-. 054782 


. 724134 

I 12697 

-.031716 1 

.9 

. (167923 

-3962 

-1. 56730 

-.165447 

061513 

-. 184643 

. 812714 

127257 

-.034084 I 

1.0 

.064712 

-4518 

-1.74432 

-. 142960 

-. 066176 

niH 

. 875829 

-. 138772 

IjgEwijBi 

1.1 

.043246 

-6321 


-. 127027 

-.063518 

-. 133012 

.908173 

-. 148055 

-.030590 , 

1.2 

.032077 

-0274 


-.090426 

-.008301 

-.081430 

. 902065 

147554 

-.034981 

1.3 

-.012565 

-8570 

-X. 75064 


-. 064585 

— . 026512 

. 856022 

-. 145166 

-. 032478 | 

1.4 

-.022802 

-6778 



-.056328 


.744686 

-. 125729 

-.025352 . 

1.5 

026023 

-5079 

-1. 25768 

.058562 


HTmI 

. 586752 

Tv^B 

-.020274 ■ 

1.6 

-.047200 

-4748 


. 091287 


.164326 


-.072582 

-.014992 ; 

1.7 

-.075428 

-4048 


iWBiKikibh^B 

-.014737 

. 213867 

. 196931 

-. 039141 

-. 008194 ! 

1.8 

071623 

-1908 

ilifiE fiWTiB 



.252820 

-.022338 

Bmh 

-.001334 , 

1.9 

-.078074 

-316 

.28440 

. 196353 


.285725 

-. 250374 

. 031215 

lIBL'^kM 

2.0 

-.082513 

1018 


Mi* 


.288821 

-. 476502 

.067093 

.011195 

2.1 

-.088003 

2172 

1. 17552 

1119 ■’! 

.051697 

. 284010 

-. 679353 

[Tn 

. 010361 

2.2 


2835 

1. 40 624 


. 065082 

. 280136 

-.854892 

.127658 

.021110 

2.3 

087504 

5070 



.076895 

. 271918 


.156658 

. 025542 ! 

2.4 

-.030227 

9801 



.085314 

. 216138 

-1. 124562 

. 183397 

. 02673 1 1 

2.5 

-.024950 



.121859 

,088583 

.125260 



I 

2.6 

-.026422 


1. 78224 

.051241 


.021081 



| 

2.7 



1. 62424 

.006889 




l 


2.8 

-.019651 


■hi 


. 066452 



1 


2.9 

-.0058G4 


1. 22608 

-.072772 


-. 147119 


1 

i 

3.0 



.99856 

-.080953 

.037762 

-. 152347 




3 1 

.005620 


.75208 


.023846 

156342 




3.2 

. 007958 


.34760 

-.083537 

,012060 

-. 136064 




3,3 

. 010611 


.29072 

-.085259 


136890 




3.4 

, 009145 


. 01264 







3.5 




062437 

-. 012536 

-.087786 
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Method B, — The principal computations illustrating 
method B are presented in table V. Again, some of the 
intermediate steps outlined in method B are simple com- 
putations and are therefore omitted. In table V(a) the com- 
putation demonstrating the determination of C„ ta , C ,„ t jf and 
C,„ s by step (4) by using the relation (25) is shown. Three 
of the columns are obtained by operating on columns ©, ©, 
and © of table III with the integrating matrix ||Cj| given 
in table I. 

Table V (b) illustrates step (7), the determination of C La 
and C Ls using equation (19). Two of the columns are ob- 
tained directly from table III and the other column is 
derived by using equation (18). 

In table V (c) final results obtained with method B for 
three sets of flight data are shown. 

Method C. — The principal computations of method C are 
presented in table VI. Table VI (a) shows the computa- 
tion of Ki, Ki, K b , and K b from flight 1 data by the method 
of reference 4. The integrals in table VI (a) were computed 
by reading the film at 0.05-second intervals and using the 
integrating matrix for Af=0.05 second; this interval was 
necessary in order to obtain reasonable results for the 
method. Use of the time interval At =0.1 second did not 
produce sufficiently accurate values of K t and K e in this 
case. Table VI (b) shows the computation of the K values 
for flight 2 data. In this case a time interval of At=0.1 
second was sufficiently small to produce reliable results for 
method C. 

In table VI (c) the final results obtained with method C 
for flight 1 and flight 2 data are given along with the results 
obtained by using the approximate formulas of appendix C. 

DISCUSSION 

The three methods presented in this report are based on the 
assumptions that the aircraft has two degrees of freedom 
(vertical motion and pitch), that the motion of the aircraft 
can be adequately described by the linear differential equa- 
tions of motion with constant coefficients based on small- 
perturbation theory, that the aircraft is a rigid body with no 
flexibility, and that the major contribution to the damping 
comes from the horizontal tail. The airplane, its flight con- 
dition, and the maneuver to be analyzed must therefore fall 
within the realm of these assumptions; that is, the airplane 
should be operating under conditions in the linear range of the 
coefficients, the maneuver should be of the pull-up or push- 
down variety where little loss in airspeed occurs during the 
maneuver and where displacement angles are small, and the 
maneuver should start from a level-flight trim condition and 
should be in the Mach number range in which these assump- 
tions are valid. 


Since the choice of the method to be used depends primarily 
on the number of measurements which are available, method 
A is recommended when four basic measurements are avail- 
able, method B when three measurements are available, ami 
so forth. If, however, an accurate value of X is known in 
advance, then method B is recommended since it will give 
the most information for the least amount of work. Method 
C requires more work than method B, and the modified 
method C is not expected to be so reliable as the other methods. 

In these methods sufficient data to cover the natural period 
of the short-period oscillations of the aircraft should be used. 
For highly damped motions sufficient data should be used to 
approach the steady-state value. 

The accuracy of the results obtained from these methods is 
influenced considerably by errors in the instruments and in 
the* record reading. Instruments used should be accurate, 
calibrated both statically and dynamically, and free from; 
drift and hysteresis. Before an analysis is started the data, 
should be corrected for known instrument errors; the records 
should then be read as carefully as possible. Measured tail- 
load data should be corrected for effects of inertia. The 
accuracy of the analysis next depends on the time interval 
selected for the integrating matrix and on the amount of 
departure from the basic, assumptions. Provided the initial 
data are accurate, the ^mailer the time interval the more, 
accurate the results. If at all possible, therefore, time inter- 
vals of Af=0.05 second or At =0.025 second should be used 
for an accurate analysis. With the introduction of oscillo- 
graph timers, which record timing marks every 0.01 second, 
and the use of IBM facilities to process the data, such timing 
intervals are feasible. The differences between the values 
shown for different flights in tables IV (e), V (c), and VI (ej, 
are believed to represent the scatter caused by effects of 
flexibility, minor nonlinearities, instrument errors, record- 
reading errors, changes in airspeed during the maneuvers, and. 
other items which essentially depart from the basic assump- 
tions. j 

\ 

As may be seen from a comparison of tables IV (a) und 
IV (d), the inclusion of the and C Lk terms in the force’ 
equation for method A has little or no effect on C L<1 but has a' 
considerable effect on If the C Lii and 0* terms arc 
retained in the force equation in the development of equal ion 
(8), the form of the equation remains the same but the K. 
values now include C L& and terms. These terms were 
found to have a negligible effect on the K values and their 
inclusion made the equations too unwieldy to handle. For 
the sake of completeness, the K values including the and 
C La terms are given in appendix E as K h K>, Ju, K u ' 
and K 6 . 



METHODS FOR DETERMINING AIRPLANE LONGITUDINAL-STABILITY DERIVATIVES FROM FLIGHT DATA 

TABLE V.— COMPUTATIONS ILLUSTRATING METHOD B 

(a) Determination of C m a , C m j, and C m , 


t 

f Aa dt 
Jo 

So A£ dt 

So AS dt 


t 

So Aadt 

J.'“* 

So AS dt 

jSc® . 

0 

0 

0 

0 

0 

1.8 

-0. 065668 

-0. 104508 

0. 046582 

0. 015219 


0 

0 

. 000182 

0 

1.9 

-. 064612 

-. 077530 

. 039040 

. 017840 

. 2 

000018 

-.001164 

. 003154 

-. 001800 

2.0 

-. 061805 

-.048742 

. 030954 

. 018857 

. 3 

-. 000295 

-.007697 

. 009716 

-.006729 

2.1 

057382 

-. 020009 

. 022443 

. 019444 

. 4 

— . 001232 

-. 022265 

. 017219 

-.011385 

2.2 

-. 051530 

. 008260 

.013272 

. 019718 

. 5 

-. 003153 

041702 

. 024501 

-. 013576 

2.3 

-. 044418 

. 035893 


. 019718 

.6 

-.006200 

-.063115 

.031585 

-. 015062 

2.4 

-. 036264 

.060589 

-. 002693 

. 016979 

.7 

-. 010314 

-.085043 

. 038487 

-.015219 

2.5 

-. 027526 

. 077951 

005296 

. 011072 

.8 

-. 015379 

-. 105918 

. 045288 

-. 014867 

2.6 

-. 018770 

. 084949 


. 004656 

.9 

-. 021209 

-. 125274 

.052139 

-.014123 

2.7 

010553 

. 083685 

010526 

. 000626 

1.0 

-. 027613 

-. 142554 

. 058791 

-. 012989 

2.8 

-. 003317 

. 076269 


-. 003590 

1. 1 

-. 034367 

-. 157253 

. 064103 

-. 011541 

2.9 

. 002640 

. 063493 

-. 014122 

-.006612 

1.2 

— . 041237 

-. 168003 

. 067784 

-. 008216 

3.0 

. 007146 

. 048510 

— . 014312 

— . 007355 

1.3 

047910 

-. 173429 

.068473 

-.004889 

3.1 

.010221 

. 033065 


-. 008059 

1.4 

-. 053989 

-.171125 

.066415 

. 001330 

3.2 

. 012003 

. 018620 

-. 013267 

-. 007596 

1. 5 

— . 059045 

161157 

.064119 

. 005321 

3.3 

. 012087 

. 005148 

-. 012304 

-. 007746 

1.6 

— . 062767 

-. 146702 

. 060601 

. 008294 

3.4 

. 012428 

-. 006739 

-. 011282 

-. 000280 

1.7 

-. 065018 

— . 127892 

. 054201 

. 012089 

3.5 

. 011447 

-. 016172 

-. 010294 

-.000869 

aSc » 

f Aa dt+Cmj 

* AE dt+ Cm, 


[ 


0. 0448788 
0. 0743828 
-0. 0370442 


0.0743828 -0.0370442-1 (C m \ ( -0.0072099 

0. 2959167 -0. 1121529 ) Cm; ( = ' 0. 0059555 


-0.1121529 0. 0471095 J (c, 


n, 


-0. 0006984 


Cm^-0. 624 
= -0.149 

Cwtj — — 0. 861 




(b) Determination of Cl„ and Ct, 


t 

Aa 

(table III, 
column ®) 

AS 

(table III, 
column ®) 

A\p 

t 

Aa 

(table in, 
column ®) 

AS 

(table III, 
column ®) 

A\p 

0 

0 

0 

0 

1.8 

0. 001884 

-0. 071623 

-0. 030434 

.1 

0 

. 009703 

0 

1.9 

. 019360 

-. 078074 

. 064502 

.2 

-.000922 

.055812 

. 012909 

2.0 

. 036467 

-. 082513 

. 192411 

.3 

-.005175 

.072880 

028063 

2.1 

. 051697 

-.088063 

. 322604 

. 4 

-.013921 

.074625 

-.091424 

2.2 

■THTP] 

-.095707 

. 407066 

.5 

-. 024871 

.071414 

-. 168381 

2.3 


-. 087504 

. 484441 

.6 

-.035932 

. 070698 

-. 253435 

2.4 

. 085314 

-. 030227 

.586746 

. 7 

-. 046216 

. 067923 

325120 

2.6 

.088583 

-. 024956 

. 530203 

.8 

-. 054782 

.068691 

-. 395488 

2.6 


-. 026422 

. 514713 

.9 

-. 001513 

. 067923 

442312 

2.7 

. 077804 

-. 026702 

. 472812 

1.0 

-.066176 

. 064712 

-. 495065 

2.8 

. 066452 

-. 019651 

.431104 

1. 1 

-. 0G8518 

.043248 

-.529718 

2.9 

. 052350 

-. 005864 

. 364215 

1.2 

-. 068301 

. 032077 

-. 501830 

3.0 

. 037762 

. 000908 

. 298661 

1.3 

-. 064585 

-.012565 

-. 505261 

3.1 

.023846 

. 005620 

. 227539 

1.4 

-. 056328 

-. 0228G2 

-. 425090 

3.2 

■L'lLhvJ 

. 007958 

. 109150 

1. 5 

-. 044150 

-. 026023 

-.372095 

3.3 

.001887 

. 010611 

. 092732 

1.0 

-. 030073 

-. 047260 

-. 266469 

3.4 

006638 

. 009145 

. 010140 

1.7 

014737 

075428 

-. 146952 

3.5 

-. 012536 

. 010175 

-. 049418 


Cc) Final results from three sets of flight data using method B with X-0.5 



Flight 1 

Probable 
error for 
flight 1 

Flight 2 

Flight 3 



7.09 

0.113 

6.68 

6.78 

CL,} 

0. 450 

0.105 

0.402 

0.420 

Cl; 

0. 062 

0.008 

0.057 

0.053 

<&* 

0.031 

0.004 

0.028 

0.026 

Cm a — — — - 

-0. 624 

0. 020 

-0. 670 

-0. 698 

Cmi 

-0. 149 

0.019 

—0. 136 

-0. 126 
-0.063 

Cm j 

-0. 075 

0.010 

-0.068 

Omg 

(OCi/ba), 

acy&s 

-0. 861 
4. 21 
2.13 

0. 063 

-0.813 

3.79 

1.87 

-0. 892 
3. 52 
1.90 


A^ = Cr„Aa+ Ci f AS 


r 0.087761 

-0. 053345 “I 

cu 

0. 598132 

L -ft. 053345 

0. 100964 J 


-0. 332334 ‘ 


Di„=7.09 

Cij=0.450 
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TABLE VI— COMPUTATIONS ILLUSTRATING METHOD C 


(a) Determination ot Si, St, St, and St from flight 1 data 


t 


X ? dt 

-J^ AS it 

- f P A Sdrdt 
jo Jo 

-A 9 

t 

A 8 

j; AS it 

-j^AS it 

- ^ f r A5 dr dt 

Ju Jo 

-Atf 

0 

0 

0 

0 

0 

0 

1.2 

-0. 133601 

-0. 0626S8 

-0. 06933 

-0. 03863 

0.090429 

.1 


0 


0 


1.3 



07124 


M W 

.2 



-.00333 


.019808 

1.4 

fllBEKS® 

-.090671 


-.05274 

-.014610 

. 3 



00933 


■tu.Ii-.lMI 

1.5 

-. 138768 

-. 104769 

-.06720 

-. 05958 

-.058562 ] 

.4 


~. 001241 

-.01734 


.125305 

1.6 

-. 131333 


-.00336 

-.06613 

— . 091287 ! 

.5 

029159 


-.02473 


. 149419 

1.7 

-. 120101 

E 

-.05716 

07218 

-.133055 

.6 

-.044962 


03205 


. 165781 

1.8 

-. 104585 

-. 142170 

-. 04975 

-.07752 


.7 

061712 

-.012469 

-.03921 

-.01008 

MmkiIMI 

1.9 

-. 086708 

-. 151759 

-. 04230 

-.08213 

-.190353 ! 

.8 

-. 078140 

019465 

-. 04024 

-.01495 

. 163628 

2.0 

-.066056 

1 r'H£SF§]'t ,Ll J,':, 1 

-. 03387 

-.08594 

-.207519 ! 

.9 

-.094122 


-.05317 

-.01992 

. 155447 

2.1 

-. 045823 

165066 

02518 

08890 



-. 109042 

-.038253 

-. 05994 

-.02558 

^HE2£IM 

2.2 

-.024343 

168577 

-.01588 



1.1 

-.122585 

-.049848 

-.06560 

-.03188 

PPInM 

2.3 

mHI i • 1 ►’>. v. i.fl 


-.00649 


-.217022 ‘ 





• 


2.4 

.018660 

-. 109008 

-.00126 

-.09211 

-.186880 


rt rt rt r r 

S,AS+S 2 I Afldt-A'. AS it—Ka ASdrdt=-AB 

Jo Jo Jo Jo 


■o. 183701353 
0. 141082469 
0.094470059 
.0.081352270 


0. 141082469 
0. 247220311 
0.073939691 
0. 130293337 


0. 094470059 
0. 073939691 
0. 049289973 
0.042731928 


0. 0813522761 / A', \ / -0. 004034232\ 

0. 136293337 ) A'j (_ ) 0. 229008158 f 
0. 042731928 -0. 008849134 ( 

0. 075308217 J \ Kt' (. 0. 11 9723023 J 


I i 

I . 

I 

I 


Si =4.14 


Kt= 9.55 
Ki=-9.77 
St= -14.62 


Cl)) Determination of Si, St, Ks, and Ki from flight 2 data 


t 

AS 

J^A9<« 

-J^ASdt 

- f f T A5 ir it 

Jo Jo 

-Ad 

t 

Ad 


- JJaS dt 

- P f r A6 dr dt 

Jo Jo 

-AS 

j 

0 

0 

0 

0 

0 

0 

1.3 

-0. 148557 

-0. 081124 

-0.075188 

-0.048259 


.1 


0 

.000094 

.000025 

0 

1.4 

-.146021 

iOiMM 


-.055559 

-.053825 

.2 



-. 002278 

-. 000063 

.013779 

1.5 


-. 110267 

-.002688 

-.002211 


.3 

-.004005 


-. 008715 

-.000599 

.065021 

1.6 

-. 122660 

-. 123344 

-.054155 

-.068062 

-.182144 

.4 

-. 013776 


-.010785 

-. 001861 

. 126160 

1.7 


-. 134634 

-.045741 


-.221750 ' 

.5 

-.028391 

-.003117 

-.025055 

003952 

. 161900 

1.8 

l - , ISM i 


-.037416 

-. 077213 

-. 230844 ! 

.6 

-.045744 

-.006811 

-.033384 

-. 00G874 

.183005 

1.9 

-. 054902 

-. 150417 

-.028839 

-.080528 

■u '/i'. iy 

.7 

-. QG45QI 

-. 012316 

-. 041523 

-.010622 

. 191186 

2.0 

031272 

-. 154720 


-. 082973 

— . 232524 1 

.8 

0837G2 

-.019738 

-.049313 

-.015167 

. 190325 

2.1 


-. 150692 

-.010793 

-. 031517 

-.238218 

.9 

mmwmm 

-.029039 

-. 056182 

-.020438 

. 174393 

2.2 

.014414 

■ m M M 


-.035120 

-.221759 

1.0 

-. 118566 


-.063462 

-.026417 

. 153294 

2.3 

.035352 

-. 153867 


084852 

-.153339 

1.1 

132553 


-.071765 

-. 033208 

.125305 

2.4 

.052386 

-. 149447 

. 012185 

-.083915 

— . 148088 

1.2 

-. 143128 

-.060497 

-.076422 

040647 

.083100 

2.5 

.065430 

143521 

.018077 

082370 

-.114540 


S 1 A 6 +S 1 f *A9 dl-Ki PaS it-Ki f r A5 (ir it=-AS 

Jo Jo Jo Jo 


r0.18385708 

0.09810197 

0.09180851 

0.05739522“! 

fA',1 

-0.02525329 1 

0.09819197 

0.25436118 

0.05665420 

0.14095293 


0.28335143 [ 

0.09180851 

0.05665420 

0.04725485 

0.03311576 

\K S \~ 

-0.01717559 

Lo.05739522 

0.14095293 

0.03311576 

0.07828319J 

UJ 

0. 14957619 J 


Ai=4.19 


A. = 10.33 
Ks= - 10.01 
St— 15.53 


(c) Final results using method C with X=0.5 



Flight 1 

Flight 2 

Accurate 

values 

Approximate 

values 

Probable 

error 

Accurate 

values 

Approximate 

values 

Cn. — 

7.21- 

^ 'l 

0.106 

7.59 

7.34 

Cl„ 

0.374 

0.371 

0.013 

0.304 

0. 391 

Cl-, 

0.066 



0.067 

0.009 

Ct-i 

0.033 

0.034 


0.034 

0.034 

Cm a ..... 

-0. 624 

-0.627 

0.073 

-0.706 

-0. 710 

Cm' — « — 

-0.158 


0. 012 

-0. 161 

-0. 164 



-0. 079 

-0.080 

■ H 

-0.081 

-0.082 

Cmj — 


-0.887 


-0.941 

-0.935 


4.46 

4.54 


4.49 

4.58 

iCLfiS 

1.75 

1.73 


1.84 

1.83 
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In the actual computation it is recommended that the 
simultaneous equations formed by the least-squares procedure 
be solved directly by the elimination of the variables or by 
Crout’s method. (See ref. 6.) The use of a least-squares 
method permits the calculation of a probable error, which is 
an indication of the fit of the data. The expression used in 
computing the probable error is * 

P.E. = 0.6745^|^Y 

where B u is the main diagonal term of the inverted matrix 
of the coefficients, E is the difference between the computed 
and measured value of the variable, N is the number of cases 
considered in the least-squares procedure, and k is the num- 
ber of unknowns determined by the least-squares procedure. 
A probable-error analysis was made of all the results using 
flight 1 data arid these results are given in tables II (b), 
IV (e), V (c), and VI (c). This probable-error analysis 
indicates that all the derivatives determined by method A 
with the exception of C Ls appear to be more accurate than 
the derivatives determined by method B or C; it also 
indicates that the derivatives determined by methods B and 
C appear to be of the same order of accuracy. 

When the computed stability derivatives are substituted 
back into the equations of motion, the method that uses the 
most measurements and has the fewest restraints imposed 
on it would be expected to produce the most accurate results 
and give the best fit to the original data. This might not 
be the case, as illustrated in figure 2 which compares the fit 
' I of the measured data with the computed data for the three 


methods presented. The results for method A are more- 
accurate for the data herein than the results for method C, 
but the fit of the incremental-pitching-velocity curve for 
method C is as good as, if not better than, the fit for method 
A or method B. It appears in general that the more coeffi- 
cients determined from a single time history the better will 
be the fit of the data but the less accurate will be the coeffi- 
cients determined. The fit of the data is interesting since 
the three methods presented are essentially curve-fitting 
processes in which the longitudinal equations of -motion are 
used to fit the flight data. A good fit indicates that the 
equations of motion and assumptions used adequately fit the 
data and the coefficients determined, if inserted in the 
equations of motion, will reproduce the motions of the air- 
craft. 

In figure 2 the incremental tail load shown for method B 
was computed by using the stability derivatives determined 
from the time histories of incremental load factor and 
pitching velocity. In method C the incremental load factor 
and tail load presented were computed by using the deriva- 
tives determined from the pitching velocity. These time 
histories indicate how well the derivatives determined on 
the basis of the measurements recorded by one set of instru- 
ments will predict the measurements recorded by a .different 
set of instruments. In the case of method C the agreement 
is good; in the case of method B it appears that a more 
realistic value of X than 0.5 should be used. Method B is 
more sensitive to X than method C is. 

Although not presented, the derivatives determined from 
flight 2 by methods A and C were used to predict the motions 
of the aircraft for flight 1. A comparison was then made 



Figure 2. — Measured and computed flight 1 time histories of incremental elevator displacement, pitching velocity, load factor, 

and tail load showing the fit of the data. 
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with the actual flight 1 motions and it was found that the 
predicted motions and the actual motions were in good 
agreement. These results verify the validity of the method 
outlined herein as applied to the example airplane. 

A possibility for a further generalized method which would 
include damping effects of whig and fuselage and therefore 
would make the method applicable to the case of sw r ept-wing 
airplanes may be realized by combining features of method C 
with method A in the following manner. Equation (A20) 
may be written in the form 

Now C m g t and C mu may be evaluated by a least-squares 

procedure, provided an accurate value of C Li f g available or 
can be determined. Examination of the probable errors for 
C Ls given in tables IV (e) and VI (c) indicate that, in the 
case of the medium jet bomber used in the calculations 
herein, the more accurate value of C Li is determined by 
method C by using equations (C12) and (C4). It is believed 
this will generally be the case for the derivative C Ly It 
might also be noted that tliis value of ft. will provide more 
rapid convergence of the iterative procedure of method A. 
The usual assumption is made that the contribution of 
wing and fuselage to C m& is negligible. Then C ttti can be 
computed through the use of the value of C m& computed by 
the above procedure and equation (A22). 

Possibilities for further investigation are to expand the 
method to include flexibility effects and the effects of higher- 
order derivatives and to extend the method to the case 
where the initial conditions are known but are not neces- 
sarily zero; that is, the maneuvers do not start from level- 
flight trim conditions. The methods could also be extended 
to other configurations such as canard aircraft and tailless 
aircraft, and perhaps a similar analysis could be made of 


the lateral motion of an aircraft to determine the lateral 
derivatives. 

CONCLUDING REMARKS 


An analysis of longitudinal-stability derivatives by three 
separate methods has been presented and applied to flight 
data. Method A, the most general method, requires four 
measurements in time-history form and permits computation 
of all the longitudinal-stability derivatives; it also requires 
the most computing time and gives the most accurate 
answers. Method B, which requires tliree measurements in 
time-history form and one supplemental assumption, namely 
C m . 

-s r -“= Constant (where C m& and are the pitehing- 

'-'mg 


moment derivatives), gives the most information for tin* 
least amount of work and gives results which are in good 
agreement with those of method A. Method C requires ■ 
two measurements in time-liistory form and two 

plementary assumptions, namely 


sup- 

~A~= Constant and 


' mg 


X 

C mt =~=-C Ls (where C ms and C Ls are the elevator-effectiveness 
c 

derivatives, x t is the tail length, and c is the mean aero- 
dynamic chord). 

The results obtained for the methods presented depend in 
a large measure on accurate instrument measurements and 
require considerable computation to yield adequate engi- 
neering answers. Since, however, the present trend is toward 
increased instrument accuracy and expanded facilities for 
machine computation, this directon appears to be the one 
in which flight-data analysis should proceed. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 15, 1952. 


APPENDIX A 


DETERMINATION OF LONGITUDINAL-STABILITY DERIVATIVES BY USING FOUR MEASUREMENTS 

IN TIME-HISTORY FORM 


EQUATIONS OF MOTION 


The equations of motion for small vertical-plane dis 
turbances may be stated (see fig. 1 for definition) as 


W 


?S4 “ + (^r), 

bC bC 

-ftp VtqS t c t A5+-^ 1 mqSiXtAS 


bC Ll n 
vtiSAS 

(Al) 


(A2) 


where 


(A3) 


These equations are for a rigid body and are based on the 
usual assumptions of linearity, small angles, and no loss in 
airspeed during the maneuver. It should be noted that the 
variables are all in incremental form measured from a steady- 
flight trim condition. 

Substituting equation (A3) into equation (Al) results in 


IT 

S 




bC Lt o 

-fig- vtqS t AS 


Since 


a V • 
A n = — 7 

g 

equation (A4) can be expressed as 
W 


where 


An — Cz Ja Aa-{-Ci, a 6i-rC L ^- J rCL S A.8 


+ 

(A4) 

(A5) 

(A6) 

(A7a) 


n — / dCV x S t x t be 

Cl *~ \ba), v ‘SVba' ' (A7b) 


a 




(bC L \ 


S V^ t 


(A7c) 


C Ll = 


*C Lt St 


bS n ‘ s (A7d) 

Substituting equation (A3) into equation (A2) gives 


.zi-iS 


s t c, . ?>C Lt s, x, (A . 
' - - r)t s ° ' a* ~ v“.8) 


b8 


bS 


S c 


or 




AS 


S t x t J>Cm t S,c,\ 

V bS v ‘ S c + bS Vt S c) 


+ 

(A9) 


which can be expressed as 
Id 


where 


^^~^”* a Aa!+ C m6c a J rC m ' l d-\-C mi A8 (AlO) 


r - §1 x * be 

\baJ t V ‘ S cVba 

f bCl ^ S t x t 2 1 


'me 


\ ba J, ' S cV 


n ■ _^ G h St x, , St C, 

bS Vt S b$ 


yvt 
S t . 

Vt ~sT 


■ (Al lb') 
(Allc) 
(Alld) 


From an examination of equations (All) and (A7), the 
following relations are seen to exist: 


xt n 

m *~C 


n _£< rt 

{ ~’ m 6 Q ■ 


(A 12a) 
(Al2b) 




(Al2c) 

Equations ,.(A6) and (AlO) are linearly dependent in the 
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present form and must be suitably altered to be put into a 
computational form. As is well-known, C L6 and C L& are 
small, and initially for computational purposes the force 
equation can be expressed as 

^An=C La Aa+C L £8 (Al3) 

The derivatives C L6 and C L& can be determined by means of 
equations (Al2b) and (Al2a), respectively, after G mi and 
Cm A are determined. 

Equations (A13) and (AlO) are now in the identical form 


of the equations of motion developed in reference 7. 
figure 1 the following relation is seen to exist: 

From 

A6=Aa-\-Ay 

(A 14) 

As demonstrated in reference 4, equations (AlO), (Al3), 
and (A14) may be solved simultaneously to obtain the rela- 

tion 


6 4- K\6 -j- K 2 A6— K 5 A8 + .Kg f Ab dt 

(A15) 

which mayHbe expresseed in integral form as 


9+KiA6 + K^‘ o A6 dt=K^ Q Abdt+K^‘J^Abd 

r dt 

where 



(A 16a) 

k,~ «f («..+«,, ;■*-) 

(A 16b) 

„ qi>C / p qS p p \ 

As— I \ 

(Al6c) 

K > =1 TmV (A.C..-GAJ 

(Al6d) 


By using the matrix method of reference 4, K u K 2 , K s , 
and K f , may be evaluated from the time-history measure- 
ments of pitching velocity and elevator angle. 


METHOD OF SEPARATING C’ m ^ AND C m& 

This method of separating C mi) and C m& applies only to 
conventional aircraft configurations equipped with a hori- 
zontal tail surface located to the rear of the, wing so that the 
major contribution to damping in pitch is due to the hori- 
zontal tail. 

In order to separate (J mi and C mu , the tail-load equation is 
developed into a form suitable for computing C mi separately 
in the following manner: ■ * 

The incremental tail load is given by 
/5>p \ bCr 

AL c =(~l)' Vt qS t Aa t + ^ Vt S t qA5 (A17) 


Substituting equation (A3) into equation (A17) gives 

^ [ A “ 0 -!)-.# * T v=]+ 


bb 


ti t S t qA5 


(A 18) 


or 


(Al'J) 


AL t _/bCr\ 8 t /. be\. fbC L \ S t X t be . 

qS~\da S \ ba) Aa \ba S V ba a 

/b(\\ S 1 x,_ + St .. 

VdaA 17 * s V-yfrt bb 1,1 S M 

which can be expressed as 

(^ t C mi -C m S\Aa+^ (C m .a+C mi 6) + r, t AS 


From equation (A 14) it can be seen that 

C n& a-\-C m fi=(C mit -\- a+f mj 7 

but from equation (A 16a) 

Cm - + Cm ^Wrn ( ’ La ~^Sc Kl=Kl,) 


(A2U) 

(A21) | 

I 

j 

(A22j ; 


Therefore, 

C mk a-{-C m gd—Kioa-{- C m gi (A23) 

Substituting equations (A22) and (A23) into equation (A20) 
gives 


~ (Kma+ Cmgi'j + C^Ab 


(A24) 


which can be expressed as 


A L ( c V 


qS x f KiaAa a :,^ I0 “ 


or 


where 


and 


=’C mi jjyY— Ip i (v'i? t +l) AaJ 
Afl = CmgAtp 


Am=~J‘ -|£ KwAa — ~ Kna-C H Ab 


(A2.<5) 

(A26) 

(A27) 


A*>=-^ T — ji (\q<+l)Aa=^|j|r An— l) AaJ 


fA28i 


From equations (A14) and (A5) the following relations arc 
self-evident: 

a=-9 — An 

Aa = A0 — yj‘ An dt (A2 r J> 

The relations needed to determine the longitudinal-stalnlh v 
derivatives from the flight measurements have now been 
developed from the equations of motion; it remains to express' 
the pertinent relations in matrix notation. 

MATRIX FORM OF THE EQUATIONS j 

A powerful tool for data analysis is provided by matrix] 
. methods since tabulated time histories may be conveniently; 
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carried in the equations of motion. In the matrix solutions 
using data, it is well known that numerical differentiation is 
inherently more inaccurate than the corresponding integra- 
tion process. For this reason, whenever necessary, the 
differential equations are expressed in integral form. The 
' first step in matrix solutions is to tabulate the recorded 
1 values of the basic variables at a number of points t 0 , t u t 2> 
t a , . . . along a given time history as in table III, the interval 
. of time used in most of the computations in the report being 
At— 0.1 second. These tabulations then become the various 
column matrices Ad { , A L, v An t , and d t . In certain cases 
■ smaller time intervals must be used to get reliable results. 
Another means of getting more accuracy is to use integrating 
matrices based on cubic or quartic curves faired through the 
data in place of the parabolic curves. 

The four basic measurements used in the development 
herein are incremental load factor, pitching velocity, tail 
load, and elevator angle. By use of equation (A29), the 
time history of incremental angle of attack is computed. 
Equation (A13) may be expressed as 


Aao 

AS 0 



'An 0 ' 

Aai 

A5i 



Ani 

A«2 

A5 2 

Kj 

w 

~g.s H 

Arii 

> 

Aa» 

A5 S 





or 

l|BH {%;!-! lte,! (A3I) 


From equation (A22) C m& is obtained as 

O m -=K 10 -G mi (A36) 

Equation (AlO) is now expressed in integral form as 

-^=6= C mot J* Aadt-f AS dt (A37) 

which can be rewritten as 

-^=6— C m6l Aa— C m ^ A0 — C?m 0 J* o Aadt+C ms A5 dt (A38) 
Now if 

A c ~~qS% ^ ~ Aot— ^ (A39) 

then a time history of Acr can be obtained and equation 
(A38) can be put in the form 


j °Aa dt 

j‘°AS dt 


r *s 

A (To 

j\a dt 

j h A 8 dt 


Aci 

j Aa dt 

J (2 A5 dt 


Acr 2 l 

J" n Aadt 

j‘ n AS dt 


A er» 

- 


lljD| i {air H A * l} (A41) 


Applying least squares, which in matrix notation involves 
premultiplication of matrix B by its transpose B ' , to equation 
(A31) yields 


(A32) 


Applying least squares to equation (A41) gives 
[D'D] j^”jj=={.D'A<rj} 
and the solution is 


(A42) 


for which the solution is 

^ = {B'An t } (A33) 

By the method of reference 4, compute K h K h K- 0 , and K a 
from the time histories of the pitching velocity and elevator 
angle. The value of K w can be obtained from equation 
(A22). Time histories of the derived Ajx and A <p functions 
can now be computed by using equations (A27) and (A28) 
since the value of C L . has been computed from equation (A33). 
Equation (A26) becomes 

; { A^i } = Cmj { Avi } (A34) 

Applying least squares to equation (A34) results in 

t ’ 

X/ (Apn A ipi) 

(A35) 

X) (A Vi ) 2 
0 


Ao t } (A43) 

In order to include the effects of the C L g and C L& terms 
initially omitted in the force equation, equation (A6) . is 
rewritten as 

A\[/= Cz, a Ac-y- Cl s AS (A44) 

where 

W 

A\l/——gAn—C L $—C L& a (A45) 

Method A may now be iterated to obtain more refined 
values of the derivatives. The values of C Li and C Lix 
determined by equations (A35), (A36), (Al2a), and (Al2b) 
are inserted into equation (A45), and a time history of A\js is 
computed. New values of C La and C L . are computed from 

[B'B]^ C c l ^ = {B'A^} (A 4 6) 
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or 


^=[B'B]^{B'A^} 

If these values of C La and C Ls are used, a new value of K l0 
and a new time history of An can be computed, which, if 
inserted into equations (A34) and (A36), yield new values 
of C mt j and C mit . The derivatives C L& and C H are again 
determined from equations (Al2a) and (Al2b) and a new 
time history of A\p is computed by using equation (A45); 
refined C La and Ol s derivatives are found from equation (A46). 
The process converges rapidly. After it has converged, 
compute C ma and C m& from equation (A43). Thus far C La , 
0 L6 , G hit , Q l j, C na , C mi , C m& , and C ms have been determined. 


Then 


\ba )t 
(Allc) as 


may be determined by rewriting equation 
(bC L \ ScV r 

s t x t ^ t 6 

An examination of equation (Al2e) shows that 

be 1 G m& 

ba C m j 

and from equation (A7d) 

s c 

b8 r, t S t H 

All the longitudinal-stability derivatives are now determined. 


APPENDIX B 


DETERMINATION OF LONGITUDINAL-STABILITY DERIVATIVES BY USING THREE MEASUREMENTS IN TIME-HISTORY FORM 

AND ONE SUPPLEMENTAL ASSUMPTION 


The three basic measurements used in method B are 
incremental load factor, pitching velocity, and elevator 
angle. The supplemental assumption made is that C m JC mi 
is a constant, that is, 


and a time history of is computed. The values of C\ tt 
and C L{ are then computed from equation (A44) 

A$=C ' La Aa-\-C l s A5 


C, n .—\C, 


me 


(Bl) 


Equation (B4) may be expressed in matrix form as 


For a first approximation the constant is assumed to be equal 
to K (see ref. 5). 

If the definition 

A$ = Xa+lJ=(l + \)0-X-f A» (B2) 

is adopted, a time history of A| may be computed. Then 

C m $ J rC mit a'=C m g Af (B3) 

The integral form of the moment equation (A37) can then 
be written 


Aa A? S J^ Ad dt (B4) 

If time histories of incremental load factor, pitching 
velocity, and elevator angle are measured and IF, q, S, V, c , 
and I are known, then equation (B4) can be put into matrix 
form and used directly to compute C ma , C m $, and C n } . 

The derivative C m& is derived from C vli by using equation 
(Bl) and G hic and C are computed from equations (Al2a) 
and (A 12b) 


These values of C La and C L g are then inserted into equation 
(A45) 

A^—^An — CLs^—O^a 


J °A a dt 

j\ dt 

j‘°AS dt 


to 

j\a dt 

JJa$ dt 

j\ S dt 


k 

J‘\a dt 

j% dt 

j\s dt 


k - 

J*"A a dt 

j‘ n A£ dt 

j‘ a AS di 


K 


or 


Applying least squares yields 


[E'E] 


fej- 

(c ms ) 


Wc im) 


and solving for the derivatives gives 



(B7) 


(BS) 
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Equation (A44) expressed in matrix form is 


Aao A5 0 


<i 

Aai A5i 


Mi 

A aj ASi 


Mz 

, , 



A a» A5„ 




Pll 


(B9) 


(BIO) 


Applying least squares to equation (BIO) gives 

\Ol 


[B'B\ = (Bll) 


and solving for the unknowns results in 
[G 1 


= (B12) 


The derivatives ) and P - are determined from 
V5a/ ( oct 55 

equations (Allc), (Al2c), and (A7d), respectively. 

Thus the derivatives c La , C Lh 0 L& , C Li , Gm a > Gmgi G m &} 

/ 5CA 5e , <>Gl 1 . 

! ( — ) i -r—i and - ; are now determined. 

V 5a / 1 5a 55 


c, 


APPENDIX C 


DETERMINATION OF LONGITUDINAL-STABILITY DERIVATIVES BY USING TWO MEASUREMENTS IN TIME-HISTORY FORM 

AND TWO SUPPLEMENTAL ASSUMPTIONS 


By the matrix method of reference 4, time-history measure- 
ments of the pitching velocity and incremental elevator 
angle are used to compute the K values from the relation 


C m =■ 


qSc 




qS I 
mVcVm( l+\) 


C L 


KJ r 

cFm(l-f-X) 

(C 8 ) 


5 K\6 K 2 AB ■= Kf,A5 r A5 dt 

0 

(Cl) 

where 

(C2a) 

(c..+o« ;•* ) 

(C2b) 

jf . Q.&Z ( n Q.S n ft \ 

j m yU Ls o m& j 

(C2c) 


(C2d) 

If it is assumed that 


, G m ^ z=z \GmQ 

and 

(C3) 

ft x t ft 

(C4) 

these six relations in six unknowns (C La , C Li , C ma , C mv 
Omg, and O m& ) can now be solved simultaneously. The follow- 
ing relations for the variables result: 


(C5) 

where 


«=“![< K+x fL 

(C6) 

and 


ft-U+x>(*J)(i& 

P) (C7) 


, I 

m ° cFm(l + X) 



mV 

qS 



m lx. 




qSc qSc 


C L 


(C9) 


(CIO) 


Approximate formulas which give a quicker evaluation 
of the derivatives with fair accuracy were derived from 
equations (C5), (C8), (C9), and (CIO) and are 


C m 


o L ~ 


mV K, 
qS K h 


/ms ~qSc 
n 1 7 /K& Tf\ 

Vn i~l+\qSc\K 5 Al ) 


-L- x — — c 

qSc K 2 Kg m 


,-fe'O 


. I r K 1_ (K£ f KxKt 1 

'qSc L 2 H-X \kJ ^(l+xjkJ 


(Cll) 
(Cl 2) 
(Cl 3) 


(Cl 4) 


The set of approximate formulas has been found to give 
results which are usually within the accuracy of the method. 
In table VI (c), a comparison is presented between results 
computed by using the approximate relations and the more 
accurate relations. The set of approximate formulas given 
by equations (Cll) to (C14) is used in the development of 
a modified method C which is given in appendix D. 



APPENDIX D 


MODIFIED METHOD C 


For some special types of longitudinal maneuvers consider- 
able information may be determined from a single time 
history. If the elevator motion is known to be of the im- 
pulse type (a blip of short duration) but its magnitude or 
time history is unknown, then the method of appendix C 
may be modified slightly to yield some of the stability 
derivatives. The method may be used with impulse-type 
forcing functions produced by ballistic devices. If the input 
is not a pure impulse but resembles one (that is, a pulse-type 
input), then the modified method may be applied after the 
elevator motion is zero. Integrals, however, must be evalu- 
ated from the zero-time trim condition but the least-squares 
procedure is applied only to the time histories after the 
elevator motion is zero. 

Since the definite integral of an impulse is a step function 
and the integral of a step function is a ramp function, let 


and 


r 

IT 


AS dt=A 
AS dr dt—At 


(Dl) 

(D2) 


Substituting these values in the integral form of equation 
(A15) which is 

.Ki A0 ~t“ .Kj (* Ad dt — jKs (* AS dt — Kt (* F ASdrdt= — 6 
Jo Jo Jo Jo 

(D3) 


results in 

Ki Ad -(- K% ^ A6 dt — K&A — K^At~ — $ 
Equation (D4) may be expressed in matrix form as 


(D4) 


A&o 

£° A9dt 

-1 

— *0 




A0i 

Ad dt 

-1 

-tl 


1 


-0, 

AO 2 

-1 

— u 


K t 

K,A 

Ia c ,v 

► = -i 

--&* - 

a e n 

£*Mdt 

-1 

~t» 



-On 


(Do i 


Equation (Do) is then used to solve for the stability co- 
efficients K \ , K- 2 , AVI, and K C A. 

The following approximate formulas presented in appendix 
C are used to compute the stability derivatives (since only 
the ratio KJK & is used, the value of A need not be evaluated ') : 


mV Iu 

' qS Ko 


ft _ . 1 I (Kt rr 

+\qSc\K s 


) 




qSc 




Kon 

A* 


Also, it is assumed that 

As indicated previously in appendix A, 


C' m& =\C mi 


I_ (K t V | KiKq 


(1 + X) Ki 


0 


C La =-C 7 , 

“ Z t 






(Dti) 

(DT) 

(D8) 

(DO) 

(DIO) 

(Dll) 


Thus the analysis of a single time history of pitching velocity can yield considerable information if it is the response to an 
elevator impulse function; however, the elevator-effectiveness derivatives cannot be found by this method. In the ease of a 
unit impulse input (A=l), K u K 2 , K h , and A 0 would be determined directly by the method of appendix D and the equations 
of appendix C could be used to determine even the elevator-effectiveness derivatives. 
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APPENDIX E 


DEFINITIONS OF STABILITY PARAMETERS 


The stability parameters of the methods presented, not 
previously defined in the original list of symbols, can be de- 
fined as follows: 
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Matrices used in the present report are defined as follows: 


Aao 

A«i 


Act2 


11 * 11 = ■ 


A a n 


M 0 

ASi 

M 2 



IIC/H is the integrating matrix given in table I, 



£ Act dt 
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£ Act dt 
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£ At dt 
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11*11= 

£ Aa dt 

£ At dt 

£ AS dt 
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£ At dt 
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